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We propose that the energy-dependent spatial modulation of the local density of states seen by 
Hoffman et al ^| is due to the scattering interference of quasiparticles. In this paper we present 
the general theoretical basis for such an interpretation and lay out the underlying assumptions. 
As an example, we perform exact T-matrix calculation for the scattering due to a single impurity. 
The results of this calculation is used to check the assumptions, and demonstrate that Quasiparticle 
scattering interference can indeed produce patterns similar to those observed in Ref. ||]. 
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Lately the possibility that a competing charge density 
wave order exists in -B^S^CaC^Os+a has received a 
lot of attention in scanning tunnelling microscopy (STM) 
studies. For example, by applying a magnetic field Hoff- 
man et al observed a four-lattice-constant modulation 
of the local density of states (LDOS) along the (0, ±1) 
and (±1,0) crystalline directions in regions surround- 
ing the cores of superconducting vortices (l). Quite re- 
cently Howald et al reported similar LDOS modulation 
in zero magnetic field for a relatively wide range of 
bias voltage Q. These results have stimulated consid- 
erable interests in the high-T c community 0]. Recently 
Hoffman et al re-examined the zero-field LDOS modu- 
lations using high resolution Fourier-transform-scanning 
tunnelling-spectroscopyQ. They found a very important 
result - the period of the LDOS modulation changes with 
bias voltage. One of us (DHL) proposed that this phe- 
nomenon can be explained in terms of the quasiparticle 
scattering interferences (QSI). 

Scattering interference is an elementary concept in 
quantum mechanics. In a nut shell, an incoming wave is 
scattered into outgoing wave by some kind of perturba- 
tion; the interference between the two partial waves gives 
rise to a spatial modulation of the amplitude of the total 
wave. LDOS modulation due to the scattering interfer- 
ence has been long observed for normal metals |5||. In the 
context of high-T c , it has been suggested that the LDOS 
modulation due to quasiparticle interference bears inter- 
esting information on the pairing symmetry of Cooper 
pairs ||. 

The purpose of the present paper is to lay down the 
theoretical basis for the QSI interpretation of Hoffman et 
aPs data. In particular we spell out the key assumption 
which allows the deduction of the normal state Fermi 
surface and the superconducting gap dispersion from the 
observed LDOS modulations]^, |j. As a sanity check, 
we perform exact (numerical) scattering calculation for 
the case of a single impurity. The results are used to 
demonstrate the validity of our assumption and that QSI 
can indeed produce LDOS modulation patterns similar to 
those observed by STM. 



Now we briefly summarize the findings of Ref. [§. By 
examining the magnitude of the Fourier transformed 
LDOS Hoffman et al found two discernible groups of 
modulation wave vectors. The first is in the (0, ±1) 
and (±1,0) directions and the second is in the (±1,±1) 
directions. The wave vectors q of the first /second 
group disperse with bias voltage V in such a way that 
|q| decreases/increases when | | increases. In addi- 
tion, for fixed V, the |q| in the first/second group de- 
creases/increases with doping. In Ref.Q this two group 
of wavevectors are identified as the momentum transfer 
when a quasiparticle is elastically scattered across the 
normal state Fermi surface as indicated by the numbered 
arrows 1 and 2 in Fig.l. Based on this interpretation 
Hoffman et al were able to deduce the Fermi surface loca- 
tion and superconducting gap function, and found fairly 
good agreement with the results of angle-resolved pho- 
toemission experiment^]. 

Formally the idea used in Ref.[|[ fjj can be precisely 
stated in the following Greens function formalism. Let 
Go(k,o;) be the (2 x 2) single-particle Greens function in 
the superconducting state in the absence of scattering, 
and T(k + q, k; uS) be the 2x2 scattering matrix for the 
k — > k + q transition. The Fourier transform of the local 
density of states is given by 



n(q,w) 
A(cbw) 



= n (q,uj) - ^[A n (q,uj) +A 22 (q, -u) 

-A* n (-q.,u) - A| 2 (-q, -w)] 
d 2 k 

-—rGo(k + q, w)T(k + q, k; w)G (k, u) 
(2tt) 2 



(1) 



In the above the subscript "0" implies the absence of 
scattering. Given Eq. ([j]) if one makes the assumption 
(the "on-shell" assumption) that the integral in A(q, w) is 
dominated by those k's that satisfy the simultaneous pole 
equations for both Go's, one obtains the result needed 
for the QSI interpretation of Ref.[|[ 0. In particular 
one concludes that at a given bias eV = w, |n(q, w)\ is 
the largest when there is a large joint density of states 
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FIG. 1: (a) Schematic plot of banana-shaped equal-energy 
contours. Contours at higher energies are not shown. The tips 
of the bananas traces the normal state Fermi surface. The ar- 
rows indicate important elastic scattering processes discussed 
in the text, (b) The energy dependence of the modulation 
wave vector qi,2 associated with the scattering processes 1 
and 2 indicated in (a). The solid lines are naive expectations, 
and the symbols are extracted from such pictures as plot in 
Fig.2. A symbol is put down for a particular energy only 
when there is a discernable peak in the wave vector space. 
It is interesting to note that only in a narrow energy range 
peaks in both direction coexist. 



associated with scattering wavevector q. 

The normal state Fermi surface of Bi-2212 consists of 
four hole-like segments. || To get a feel of the joint den- 
sity of states let us first concentrate on the curves of 
constant energy (CCE) for the quasiparticles. (Since ex- 
perimentally LDOS modulation is only observed for en- 
ergies lower than ~ 30 meV, we shall concentrate on that 
energy range here.) Around each of the four gap nodes 
the CCE evolves from a single point at zero energy, to 
banana-shaped closed contours at higher energies. The 
size of the banana increases with the energy until its tips 
touch the Brillouin zone. When that happen the CCE 
changes from closed to open contours. This is schemati- 
cally shown in Fig.l, where different CCE's represent dif- 
ferent quasiparticle energies. Due to the large difference 
between vp and «a (vf/va ~ 20 at the gap nodes, Vp — 
the normal state Fermi velocity, and va — the derivative 
of the superconducting gap along the direction tangen- 
tial to the Fermi surface) the CCE moves with energy the 
fastest at the tips of the banana. As a result we expect 
those q's connecting the tips to have the largest joint 
density of states. This argument leads to the prediction 
that the most prominent LDOS modulation wave vectors 
should be the set of vectors connecting pairs of banana 
tips§. 

As we pointed out earlier, the above expectations are 
based on the on-shell approximation. The result of an ex- 
act scattering calculation using Eq. (|l]) can in principle 
deviate from the above expectations due to the follow- 
ing complications. 1) The "off-shell" contributions: i.e. 
those k's that do not fulfill the simultaneous pole equa- 
tions of the two Gq 's in Eq. ([!]) . These off-shell contribu- 



tions act to blur the modulation wave vectors predicted 
above. 2) The effects of scattering matrix element. Even 
when the on-shell approximation works well, the inten- 
sity of the actual LDOS modulation can depend on the 
actual scattering matrix clement, which in turn depends 
on the coherence factor. For example when the scatter- 
ing is due to a single scalar (V s ) and/or magnetic (V m ) 
impurity, the k — > k' scattering matrix element is pro- 
portional to (V m + Vs)u k u k i + {V m — V s )v k Vk' , wit h u k = 
±sgn(A fc )V(l ± e k /E k )/2,v k = Te k /E k )/2. In the 
above the upper/lower sign applies for positive/negative 
energy quasiparticle states, e k is the normal state disper- 
sion, Afe is the gap function, and E k — t\ + . For 
d-wave pairing A k changes sign in some of the scatter- 
ing processes and not in others. For example in Fig.l 
A*, changes sign in the process labelled by the numbered 
arrow 1, while maintains the same sign in the process 
labelled by the numbered arrow 2. Therefore depend- 
ing on whether the scattering process is caused by scalar 
or magnetic impurity processes 2 or 1 will be relatively 
suppressed. 

In the following we perform an exact calculation where 
the scattering is caused by a single impurity. The purpose 
of this calculation are two folds: 1) through an exact cal- 
culation we can check the on shell assumption that is cru- 
cial to interpretation of Hoffman et al's data. 2) Through 
the exact calculation we can find out whether the effect of 
scattering matrix element can make the scattering inter- 
ference patterns different from those expected from the 
joint density of state argument. This type of calculation 
(i.e. an impurity in an d-wave superconductor) is cer- 
tainly not new flllfl . However as far as we know it is the 
first that demonstrate QSI can indeed give rise to LDOS 
modulation patterns observed by STM[Q, fjj. Here is the 
summary of our results: 1) we find that among the possi- 
ble wavevectors the most pronounced LDOS modulation 
are associated scattering processes 1 and 2 (and their 
symmetry equivalence) in Fig. 1. The Fourier transforms 
of the LDOS are given in Fig.2. A comparison between 
the q — uj dispersion expected from the joint density of 
states argument and that extracted from Fourier trans- 
forming the exact LDOS modulation is given in Fig. 1(b). 

For the quasiparticle Hamiltonian, we use a model pro- 
vided by Norman et a/pO||: 



kcr 



c kn — 



k 



(A fe C fcT C_ fci + h.c). (2) 



In the above e k = J2 n =o tn Xn(k) and A k = A (cosk x - 
cosfcy)/2, where t -s = 0.1305, -0.5951, 0.1636, 
-0.0519, -0.1117, 0.0510 (eV), and A = 0.025eV. 
Moreover Xo-s{k) = 1, (cosk x +cosfc. y )/2, cos k x cos k y , 
(cos2fc 2 . + cos2k y )/2, (cos 2^;,, cos /c y + cos 2k y cos k x )/2, 
and cos 2k x cos 2k y . The Go(k, ui) and T(k + q, k;w) 
[= T(u>) for a single impurity] in Eq.(Qj) are given by, 
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For the scalar case, 
please see scalar.jpg. 



FIG. 2: Fourier amplitudes of the LDOS maps at the given 
energies as a function of momenta in the first Brillouin zone 
for (V 3 , Vm) ~ (WOmeV, 0). The numbers on the sub-graphs 
show the energy in units of meV. The hot /cold color denotes 
higher/lower strength. 



For the magnetic case, 
please see magnetic.jpg. 



FIG. 3: The same as Fig.2 for (V s ,V m ) = (0, WOmeV). 

respectively, 

Gg 1 (k, u) = (uj + iS)I - e k CT 3 - A k cri 

d 2 k 



r->) = (V s a 3 + V m I)- L - 



G (k,w),(3) 



where Cj's are the Pauli matrices. 

In order to accurately take into account the contribu- 



tion to Go from the thin bananas in Fig. 1(a), a momen- 
tum resolution 5k -C ttva/vf must be achieved. This 
in turn means that a large real space lattice has to be 
used for the actual calculation. For the results reported 
below a 400 x 400 lattice is used. In this calculation we 
did not take the self-consistent suppression of the pair- 
ing amplitude near the impurity into account. However 
we do not expect this omission to cause any significant 
change because the impurity scattering strength we used 
is much larger than the maximal pairing gap.|l2| 

In the following two kinds of impurities are studied: 
(a) non-magnetic impurity (V s — lOOmeV and V m = 0) 
and (b) magnetic impurity (V s — and V m — lOOme^). 
By using two extreme types of scatterer we can study 
the effect of matrix element on the QSI pattern. The 
impurity strength we choose is intermediate in the sense 
that we do not see quasiparticle resonance states in the 
LDOS.[|l3) It is interesting that although this type of im- 
purity does not show up in the zero bias resonance image, 
they do have important effect on the LDOS modulation. 

Fig.2 shows the Fourier amplitude of the LDOS maps 
at oj = —24 — > —3 meV and ui = 3 — > 24 meV with 
3 meV interval for a nonmagnetic. In these figures the 
intensity at q = is subtracted so that weaker features at 
other wave vectors can show up. In constructing Figure 
2 the LDOS of the central 51 x 51 plaquette in the 400 x 
400 lattice are Fourier transformed. It is interesting to 
note that in real experiment a large field of view is also 
employed in the Fourier transform in order to achieve 
high momentum resolution. 

For the negative energies in Fig.2 the Fourier peaks in 
the (±1,±1) direction is clearly visible between —3 to 
— 12 meV. Moreover as the binding energy increases the 
peaks moves away from the origin. For binding energy 
18 meV and above the Fourier peaks in the (0, ±1) and 
(±1,0) appear. Contrary to the diagonal direction spots 
these peaks move only slightly as the energy varies. Care- 
ful analysis of such movement indicates that these peaks 
move toward the origin as the binding energy increases. 
The energy dependence of the these LDOS Fourier peaks 
is shown in Figure 1(b) as symbols. The difference be- 
tween the naive expectation (solid lines) and the results 
of T-matrix calculation (symbols) is very small. This 
comparison clarifies the extent by which the qualitative 
picture of joint density of states works. In addition to 
the above peaks there are other weaker features in Fig.2. 
These features can either be identified as the higher har- 
monics of the diagonal and vertial/horizontal spots or 
can be attributed to the other scattering processes not 
shown in Fig.lQ. Upon the reversal of the bias voltage 
we find the diagonal Fourier peaks become much weaker. 
On the contrary the (0, ±1) and (±1,0) direction peaks 
seem to be insensitive to the bias reversal [|l4"|. 

The quasiparticle interference patterns for magnetic 
impurity in Fig. 3 are quite different from those shown 
in Fig.2. For example the amplitude of the (±1,±1) di- 
rection modulation for \lo\ < 12 meV is much weaker. 
The fact that the (±1,±1) direction peaks are weaker 
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than the (±1, 0) and (0, ±1) direction peaks is consistent 
with the coherence factor effect discussed previously. 

In addition to the above differences there are similar- 
ities between Figs. 2 and 3. For example at |cj| = 21,24 
meV, the patterns in these figures become rather simi- 
lar. In the same energy range the asymmetry between 
positive and negative bias becomes small. In this en- 
ergy range the horizontal component of the scattering 
wavevector associated with process 2 cease to change 
with energy. As the result the modulation wave-vectors 
in the (±1,0) and (0, ±1) directions hardly change. 

Upon closing a few remarks are in order. 1) The 
present result demonstrates that QSI due to scattering 
by scalar (non-unitary) impurities can produce LDOS 
modulation similar to that observed in recent STM 
experiment This effect is due to quantum inter- 

ference, and does not require the presence of a charge 
density wave order. 2) Whether the phenomenology of a 
very weak/fluctuating charge density wave order is con- 
sistent with the experimental observations^ ^] should be 
further investigated (|, [l5|. 3) The results of the present 
paper is obtained assuming sharply defined quasiparti- 
cles in the superconducting state |lq|. However our gen- 
eral formulation (Eq. (|l])) allows us to treat not so well 
defined quasiparticles as well - all we have to do is to 



modify the sharp pole structure of Go- We find that 
the (±1, 0), (0, ±1) direction modulation is robust against 
quasiparticle life time broadening whereas the (±1,±1) 
direction modulation is not. 4) In the present paper we 
use a single impurity as a representative of elastic quasi- 
particle scatterer. However we believe that the QSI idea 
should work under wider circumstances. 5) It is not clear 
to us what is the main source of quasiparticle scattering 
in BSCCO. For example, recently Podolsky et al argued 
that scattering by a disordered charge density wave could 
give rise to similar energy dependent (±1,0) or (0, ±1) 
modulations as well|l7j|. 6) Finally, what is the rela- 
tion between the relatively weak quasiparticle interfer- 
ence observed at zero magnetic field Q |?| and the rel- 
atively strong checkerboard LDOS modulation near su- 
perconducting vortices Q is currently unclear. 
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